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A  Class  of  Energy  Levels  for  the  Heisenberg  Linear  Chain.  II. 
Levels  near  the  Antiferromagnetic  Ground  State. 

Robert  B.  Griffiths* 


Abstract 

The  distribution  of  "unbound"  states  in  the  Bethe  formalism  for  the 
linear  chain  of  spin  £  atoms  with  a  Heisenberg  exchange  interaction  between 
nearest  neighbors  is  investigated  In  the  vicinity  of  the  antiferromagnetic 
ground  state,  utilizing  the  spin  wave  states  discovered  by  des  Cloizeaux 
and  Pearson.  An  upper  bound  is  obtained  for  the  free  energy  of  the  anti- 
ferromagnetic  chain  at  very  low  temperatures.  Plausible  arguments  are  pre¬ 
sented  to  show  that  the  "unbound"  states  actually  make  a  negligible  contri¬ 
bution  to  the  partition  function  of  an  infinite  chain,  at  least  in  the 
absence  of  a  magnetic  field. 
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I.  Introduction. 

The  Hamiltonian 

ti 

71  9  5*  *  .§  *+i 

£*■*1  *  •• 

for  a  linear  chain  of  spin  £  atoms  closed  upon  itself  (S^  is  the  spin 

operator  for  the  i'th  atom  and  J  the  "exchange  integral")  has  been  the 

subject  of  several  investigations.  The  entire  spectrum  for  a  few 

fairly  short  chains  (up  to  N*10  or  11 )  has  been  calculated  by  Orbach1, 
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Bonner  and  Fisher  ,  and  the  author  .  For  the  infinite  chain  very  little 
is  known  about  the  distribution  of  levels  lying  between  the  ferromagnetic 
(J  <  0)  ground  state  and  the  antiferromagnetic  (j  >  0)  ground  state.  The 
energy  of  the  latter  was  calculated  by  Rulthen*  using  a  formalism  due  to 
Bethe 

In  a  previous  report^,  hereafter  referred  to  as  CELI,  we  derived  the 
maximum  and  minimum  energies  as  a  function  of  the  total  spin  of  the  state, S, 
for  levels  belonging  to  the  class  C  of  "unbound  states"  in  the  Bethe 
formalism.  (For  a  brief  summary  of  Bethe' s  procedure  and  a  precise  defini¬ 
tion  of  class  C,  see  CELI,  Sec.  II.)  Recently,  des  Cloizeaux  and  Pearson^ 
have  calculated  the  energies  of  the  lowest  lying  states  in  the  antiferro¬ 
magnetic  chain  as  a  function  of  wave  vector  q.  These  "spin  wave"  states, 
as  well  as  the  ground  state,  belong  to  class  C.  The  present  report  is  an 
extension  of  des  Cloizeaux  and  Pearson's  work  to  a  detailed  study  of  the 
distribution  of  levels  in  class  C  near  the  antiferromagnetic  ground  state. 

When  we  began  the  calculation,  we  hoped  to  derive  the  low  temperature 
properties  of  the  chain  using  only  states  in  class  C.  This  category  contains 


(1) 
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a  small  fraction®  of  all  the  states  (20,^H  vs  a  total  of  2®),  but  a 
sufficient  number  to  yield  a  free  energy  proportional  to  N.  Furthermore, 
some  of  the  lowest  levels  of  the  antiferromagnetic  chain  belong  to  C. 

However,  our  calculations  lead  to  the  conclusion  that,  in  the  absence  of  a 
magnetic  field  (and  in  the  limit  of  an  infinite  chain),  the  thermal  proper-  . 
ties  are  determined  entirely  by  the  levels  in  class  B,  the  complement  of  C, 
at  all  temperatures . 

This  essentially  negative  conclusion  is  reached  via  two  routes. 

Section  VII  contains  the  direct  route:  a  very  plausible  argument  based 
upon  magnetic  susceptibility.  The  more  leisurely  route  begins  in  Section  II 
with  an  "educated  guess"  for  the  energy  of  the  chain  at  low  temperatures, 
which  is  compared  with  the  usual  spin-wave  theory.  Section  III  begins  the 
discussion  of  class  C  with  an  extension  of  dee  Cloizeaux  and  Pearson's 
results  to  states  containing  several  spin  waves.  The  free  energy  for  the 
modified  spin-wave  theory  thus  obtained  is  calculated  in  Section  IV  and 
shows  fair  (though  spurious.')  agreement  with  the  results  of  Section  II. 

An  unsatisfactory  feature  of  spin  waves — the  Introduction  of  spurious 
states — is  remedied  in  Section  V.  Several  independent  checks  Indicate  that 
the  "balls  in  pockets"  model  there  introduced  provides  a  very  good  approxi¬ 
mation  to  the  lowest  levels  in  class  C.  By  means  of  the  model,  the  contri¬ 
bution  of  states  in  class  C  to  the  total  partition  function  at  low  temperatures 
is  estimated  in  Section  VI.  The  resulting  upper  bound  on  the  free  energy 
(lower  bound  on  its  absolute  value)  is  definitely  in  excess  of  the  estimate 
in  Section  II,  but  may  still  be  useful  for  some  purposes.  Section  VIII  con¬ 
tains  a  summary  of  our  conclusions. 
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II.  Lov  Temperature  Properties  from  Finite  Chains 

For  short  chains  the  various  thermal  properties  may  be  obtained  by 
inserting  the  eigenvalues  of  the  Hamiltonian  (l)  (obtained  by  use  of  a 
high  speed  computer)  into  a  partition  function 


o 

Bonner  and  Fisher  calculated  the  energy  and  magnetic  susceptibility  for 
H»2,  3,  .  .  11;  the  results  for  N=2,  3,  .  .  .,10  were  also  obtained 

by  the  author.  ^  We  sire  much  Indebted  to  Dr.  Fisher  for  Bending  us  his 
results  for  the  case  N=ll  prior  to  publication,  and  we  have  made  use  of 
these  in  Figs.  1,  2,  and  5* 

The  calculated  thermodynamic  quantities  for  finite  chains  show  a  very- 
regular  behavior  as  a  function  of  H.  If  the  regularities  present  for 
N  s  11  persist  for  N  >  11,  the  calculations  provide:  1.  Very  good  esti¬ 
mates  for  energy,  entropy,  and  magnetic  susceptibility  at  temperatures 
above  j/k .  2.  Upper  and  lower  bounds  for  the  same  quantities  at  all 

temperatures. 

The  energy  U  of  the  linear  antiferromagnetic  chain  varies  as  T2  at 

o 

low  temperatures  according  to  the  usual  spin-wave  theory.  Let  us  assume, 
more  generally,  the  behavior: 

U--T"  (3) 

with  a  and  n  to  be  determined  by  calculations  on  finite  chains. 
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The  corresponding  entropy  ia: 

S/k.  *  W 

We  shall  henceforth  assume  that  ai  i  extensive  thermodynamic  variables 
have  been  normalized  by  dividing  them  by  the  number  of  spins  in  the  chain. 

In  Fig.  1,  energy  (measured  from  the  exact  ground  state  of  the  infinite 
chain)  is  plotted  as  a  function  of  entropy  for  chains  containing  10  and  11 
atoms.  Note  that  for  N  odd  the  ground  state  is  degenerate  and  hence  the 
entropy  is  not  zero  at  zero  temperature.  The  energy  vs.  entropy  curve  for 
N=8  (not  shown  in  Fig.  1)  lies  below  the  curve  for  N=10,  that  for  N=6  below 
that  for  N=8,  etc.  The  curve  for  N=9  lies  above  the  curve  for  N=ll,  that 
for  N=7  above  that  for  N=9,  etc.  These  regularities  suggest  that  the 
energy  vs.  entropy  curve  for  an  infinite  chain  lies  between  the  curves  for 
N=10  and  11. 

The  dotted  curve  in  Fig.  1  represents  the  function 

U/J  -  1.19  tSM"  (5) 

Curves  with  the  exponent  in  (5)  equal  to  1.8  or  2.0  which  pass  midway  between 
the  N=10  and  11  curves  at  S/k=.20  in  Fig.  1  are  definitely  low  and  high, 
respectively,  at  S/k=.30.  Let  7,  the  normalized  temperature,  be  equal  to 
kT/j.  The  estimates 

U/J  =  .V  TU 

S/k.  -=  .4-0  T  u  (6) 

follow  immediately  from  (3),  (4),  ( 5 )  and  provide  an  "educated  guess”  for 
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the  low  temperature  behavior  of  these  quantities. 

Pig.  2  shows  the  energy  as  a  function  of  temperature  for  chains  con¬ 
taining  10  and  11  spins.  The  dotted  curve  is  a  plot  of  Eq.  (6).  The 
curve  for  N=8  (not  shown)  lies  below  that  for  N=10;  the  curve  for  N*9  above 
that  for  N=ll,  etc.  This  suggests  that  the  energy  for  an  infinite  chain 
should  lie  between  the  curves  N=10  and  11. 

Eq.  (6)  shows  surprisingly  good  agreement  with  the  T2  dependence  for 
the  energy  predicted  by  spin-wave  theory.  We  do  not  regard  the  difference 
in  exponents — 2.1  instead  of  2 — as  very  significant,  since  (6)  has  been 
obtained  by  extrapolating  results  which  are  less  precise  as  the  temperature 
decreases.  3he  calculations  in  Section  VI  below  indicate  that  if  (3)  is 
the  correct  asymptotic  form  for  the  energy,  the  exponent  n  must  be  less 
than  or  equal  to  2.  Hence  it  seems  likely  that  the  T2'1  behavior  predicted 
by  (6)  is  gradually  replaced  by  a  T2  behavior  for  very  low  temperatures. 

The  following  values  for  energy,  entropy,  and  free  energy  are  then  not 
unreasonable 

U /J  .IS  7X 

S/k  »  .3 hr 

r 

f/j  «  -./*7*  (7) 

The  coefficients  have  been  chosen  to  make  the  energy  predicted  by  (6)  and 
(7)  approximately  equal  in  the  vicinity  of  7  *  *2  or  .3.  It  should  be 
eaphasized  that  (7)  represents  an  educated  guess,  and  the  values  must  be 
used  with  caution,  nevertheless,  it  appears  that  the  prediction  of  spin-wave 
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theory^  for  the  coefficient  of  t2  in  the  expression  for  the  energy  at 
low  temperatures,  n/6  =  .52,  is  definitely  too  large  by  a  factor  of 
2  or  3.  This  is  borne  out  by  Figs.  1  and  2  which  show  the  spin-wave 
predictions  for  energy  as  a  function  of  entropy  and  temperature. 
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III.  Spin-Wave  Model  for  Class  C. 

In  Be the ' formalism,  an  eigenstate  of  the  Hamiltonian  (l)  with 
total  spin  S  and  wave  number  q  is  associated  with  a  set  of  ^N-S 
integers  A-  satisfying 

0  <  Aj  H 

=  X'tT  **i  /n 

where 

M  *  Z  A;  N)  (10) 

j 

The  class  C  contains  "unbound"  states  which  correspond  to  sets  of  integers, 
no  two  of  which  differ  by  less  than  2;  all  other  states  are  in  class  B. 

(For  a  more  detailed  discussion,  see  CELI,  Section  II.)  The  ground  state** 
(S=0)  corresponds  to  the  set  of  integers 

1,  3,  5,  •  .  H-l. 

7 

Recently  des  Cloizeaux  and  Pearson  have  shown  that  the  lowest  excited 
state  (S=l)  of  wave  number  q  corresponds  to  the  set: 

1,3,5,  ...,(N-2m-l),  (N-2BH-2),  (N-2mf4), . . .,  (H-2)  for  q>0 

(11) 

2,4,6, ...,(2|m|-2),  (2|m|+l),  (2|b|+3), . . .,  (»-l)  for  q  <  0 
where  q  and  m  are  related  by  (9).  The  energy  of  the  state  (11)  minus 


(8) 

(9) 
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the  energy  of  the  ground  state  Is: 

=  ir  T  (sin  (12) 

provided  N  is  large. 

Pig.  3  illustrates  the  nature  of  the  sets  f  Aj  ^  for  N=l6.  Each  row 
corresponds  to  one  eigenstate  of  the  chain.  When  some  A  is  equal  to  an 
integer  k,  a  box  is  placed  in  the  k'th  column.  Columns  in  which  no  box 

f 

appears  contain  dots  of  "spacers."  Each  pair  of  boxes  is  separated  by  at 
least  one  spacer  for  states  belonging  to  class  C.  State  A  is  the  ground 
state.  State  B,  the  lowest  excited  state  for  m=3  (or  q=3(2n/l6) ),  may  be 
produced  from  A  by  inserting  a  spacer  in  column  11  and  displacing  the 
pattern  in  columns  11  to  15  one  step  to  the  right,  discarding  the  portion 
displaced  to  the  right  of  column  15.  An  anlogous  procedure  results  in 
state  D,  the  lowest  excited  Btate  for  m=-3- 

If  one  inserts  an  additional  spacer  in  column  7  of  state  D  ami  shifts 
the  pattern  in  columns  1  to  7  one  step  to  the  left,  the  resulting  state  E 
corresponds  to  two  "spin  waves"  with  m=-3>  The  second  spin  wave  also 
corresponds  to  ro=-3  (and  not,  for  example,  m=-4)  because  three  boxes  were 
shifted  to  the  left  when  the  spacer  was  added  in  column  7-  See  Eq.  (10). 
States  with  several  spin  waves  are  obtained  by  repeating  the  procedure  Just 
described. 

The  energy  of  a  state  containing  two  spin  waves  is  equal  to  the  sum 
of  the  energies  of  the  corresponding  single  spin-wave  states  (given  by 
(12)),  plus  a  correction  term  of  order  l/H.  An  analogous  result  bolds  for 
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a  state  containing  several  spin  waves,  provided  the  total  number  is  snail 
compared  to  N.  Consider,  for  example,  a  state  with  two  spin  waves  and 
assume  the  wave  numbers  and  are  both  negative.  Eq.  (19)  of  referenoe  7 
is  replaced  by: 

AW*  x  +  -»)-©{  '-t-'-*)- 0(^+4-*)]  03) 

The  reader  may  verify  that  des  Cloizeaux  and  Pearson's  procedure  applied 
to  the  modified  A(x)  in  (13)  yields  the  result  stated  above. 

States  produced  from  different  combinations  of  spin  waves  axe  not 
necessarily  independent.  With  reference  to  Fig.  3  suppose  that,  starting 
from  the  ground  state  we  add  one  qp  in  wave  with  m=8  (state  F)  and  another 
with  m=5.  Ifce  end  result,  state  0,  is  identical  to  state  D,  which  contains 
a  single  spin  wave  with  m=-3.  A  little  reflection  shows  that,  in  fact,  all 
states  in  class  C  may  be  produced  by  superimposing  only  spin  waves  of  posi¬ 
tive  (or,  if  one  prefers,  negative)  q.  Roughly  speaking,  we  cm  produce  a 
state  with  q  negative  by  superimposing  two  states  of  positive  q,  one  of 
which  has  wave  number  nearly  equal  to  n. 

From  this  point  on  we  shall  restrict  our  considerations  to  spin  waves 
with  positive  q.  Even  such  states  need  not  be  independent.  State  H  in 
Fig.  3,  for  example,  may  be  thought  of  aB  composed  of  either  two  spin  waves 
with  m=l  and  one  with  m=2,  or  four  spin  waves  with  m*l.  However,  if  there 
are  only  a  few  spin  waves  present,  there  is  an  approximate  agreement  between 
the  number  of  spin-wave  states  and  the  number  of  states  in  class  C,  as  we 


now  show. 
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The  number  of  states  in  class  C  with  a  given  total  spin  S,  correspond- 

g 

ing  to  ^N-S  non-zero  A  values,  is  given  by: 


n(s)  =  US+l) 


iv  +5 
1 5 


(14) 


States  with  one  or  two  spin  waves  present  are  specified  by  ^N-l  non-zero 
A  's  and  hence  correspond  to  S=l;  states  with  three  or  four  spin  waves  to 
S=2,  etc.  For  a  given  S  the  number  of  spin-wave  states  iB  thus: 


n  (s)  =>  Has)  + 


(15) 


where 


fU)* 


iN  +  y-i 


(16) 


is  the  nuiriber  of  states  containing  b  spin  waves,  assuming  the  spin  waves 
obey  Bose  statistics. 

Let  S=^pN.  One  obtains  by  use  of  Stirling's  approximation: 


lo)  Pi(S)  -  Uj  If  ■+<?((')]  +(T(Ui  N) 

l»j  « (S)  *  pN  [l-  ly  If  -hp  +■  0(pl)  ]  +  O’  (/^  m) 


(17) 


Eq..  (17)  shows  in  what  sense  n(S)  and  5(S)  are  "approximately  equal"  for 
11 


small  values  of  S. 
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IV.  Partition  Function  for  Class  C:  Spin-Wave  Model 

One  can  estimate  the  contribution  of  states  In  class  C  to  the  total 
partition  function  using  the  spin-wave  approximation  developed  in  Section 
III.  Before  doing  so  we  make  a  slight  digression  to  discuss  the  signifi¬ 
cance  of  such  a  calculation.  The  partition  function  for  the  antiferromag¬ 
netic  chain, 


Z  «  I  e_|SE" 

M 


(18) 


where  g  =  l/kT  and  the  Eq  are  the  energy  levels  of  the  chain,  may  be 
written  as 


I  *  Z8  4-  zc  (19) 

where  Zg  is  that  part  of  Z  which  comes  from  a  sum  over  states  in  class  B, 
and  Zg  that  part  coming  from  a  sum  over  states  in  class  C.  The  total 

Q 

number  of  states  in  class  C°  is  very  much  less  than  the  total  number  of 
states  for  the  chain.  On  the  other  hand,  as  noted  above,  certain  states 
in  class  C  are  lover  in  energy  than  the  corresponding  states  in  class  B, 
and  hence  it  is  not  inmediately  evident  whether  or  not  Zc  makes  a  signifi¬ 
cant  contribution  to  the  total  partition  function  Z. 

The  free  energy  (per  spin)  for  the  chain  is  given  by: 

r-  -(*(*)•'  ly  2  <»> 

Define: 


K  -  - 


(**)•'  Uj  2C 


(21) 
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Clearly, 


Fc  z  F  (22) 

or.  In  other  words,  Fc  provides  an  upper  bound  on  the  free  energy.  Under 
certain  circumstances,  Fg  might  be  equal  to  F  in  the  limit  of  large  H; 
for  example.  If  there  were  a  constant  ^  independent  of  H  such  that 

Zc/Z  >  M  >  0 

as  H  becomes  infinite .  In  the  remainder  of  this  report  we  shall  occasion¬ 
ally  use  Fg  as  a  convenient  measure  for  the  magnitude  of  V  and  refer  to 
it  loosely  as  the  "free  energy." 

Let  jus  calculate  Fg  using  the  spin-wave  approximation.  19m  problem 
is  merely  that  of  obtaining  the  free  energy  for  a  Bet  of  M/2  uncoupled 

]P 

harmonic  oscillators  with  energies  given  by  Eq.  (12)  j  hence 

,Vv 

Fc  ®  (zfl)  1 1*  !•}  [l-  J  *'«  t (23) 

Let  the  normalised  temperature  7  be  equal  to  (pj)”1.  Then  for  small  T 
the  asymptotic  value  of  (23)  is: 

P*/:r~-7V*  (2^ 

At  7  -.6,  Bq.  (24)  gives  Fg/j  ■  -.06,  vbsreas  numerical  integration  of 
(23)  yields  -.062.  At  lower  temperatures  the  asymptotic  estimate  Is  even 
better. 
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The  coefficient  of  JZ  in  (24)  it  soswwhet  lets,  In  absolute 
mgnltude,  than  the  coefficient  in  (7) •  However,  as  ve  shall  show  in 
Section  VI,  the  spin-wave  approximation  actually  leads  to  a  significant 
overestimate  of  the  contribution  of  states  In  class  C  to  the  partition 
function,  and  hence  even  the  approximate  agreement  between  (24)  and  (7)  Is 
spurious. 
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V.  "Bolls  in  Pockets"  Model  for  Class  C 

The  spin-wave  approximation  to  the  states  in  class  C,  discussed 
in  Sections  III  and  IV,  suffers  from  at  least  one  serious  defect:  the 
introduction  of  spurious  states.  It  is  true  that  if  the  number  of  spin 
waves  is  small  the  "overcounting"  is  not  large;  this  is  shown  by  Eq.  (17) • 
Nevertheless,  one  suspects  that  the  most  serious  overcounting  occurs  for 
states  in  which  the  spin  waves  all  have  low  energies,  and  these  states, 
of  course,  make  an  important  contribution  to  the  partition  function  at 
low  temperatures.  We  introduce  the  following  model  to  avoid  overcounting. 

Consider  a  potential 


V  W  =  -IxT  Sin  X  (25) 

for  x  between  0  and  n.  Let  the  interval  from  tt/2N  to  rr(l-l/2N)  be 
divided  into  N-l  equal  subintervals  or  "pockets,"  as  shown  in  Pig.  4  for 
the  case  N=l6.  Some  of  the  pockets  contain  balls,  which  correspond  to 
the  boxes  of  Section  III  and  Fig.  3.  Each  set  of  m  A  's  for  a  state  in 
class  C  corresponds  in  an  obvious  manner  to  a  distribution  of  m  balls 
in  pockets  satisfying  the  two  rules:  1.  Each  pocket  is  empty  or  contains 
one  ball.  2.  Two  adjacent  pockets  cannot  both  contain  balls. 

A  ball  in  a  pocket  centered  at  x  shall  have  a  potential  energy 
V(x),  and  we  define  the  energy  of  a  state  to  be  the  total  potential  energy 
of  the  corresponding  distribution  of  balls  in  pockets.  The  ground  state 
of  the  system  (illustrated  in  Pig.  4  for  5*16),  for  which  M/2  pockets  are 
filled,  has  the  energy 
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£,  IT  N  (26) 

if  N  is  large. 

A  state  in  class  C  characterized  by  r  =  ^N-S  A  values  has  a 
degeneracy  2S+1,  whereas  it  corresponds  to  a  single  state  in  our  model. 
However,  the  factor  2S+1  is  of  no  consequence  in  computing  the  partition 
function  in  zero  magnetic  field  (see  the  appendix  for  further  remarks) 
and  hence,  for  this  purpose,  the  model  provides  essentially  exact  counting 
of  states. 

The  separation  between  the  highest  and  lowest  energy  states  in  our 
"balls  in  pockets"  model  is  2JN,  whereas  in  class  C  (which  contains  both 
the  highest  and  the  lowest  levels  of  the  linear  chain)  it  is  I.386  JH. 

Thus  the  spectrum  of  the  model  cannot  coincide  exactly  with  the  spectrum 
of  class  C.  nevertheless,  the  model  appears  to  give  a  good  approximation 
for  the  low-lying  levels  in  class  C,  as  we  shall  now  show. 

In  the  first  place,  as  the  reader  may  verify,  the  model  correctly 
reproduces  the  energies  (see  Eq.  (12))  of  the  single  spin-wave  states 

7 

discussed  by  des  Cloizeaux  and  Pearson;  in  fact,  the  potential  (25)  was 
chosen  for  Just  this  reason.  Similarly,  for  a  state  where  only  a  few  spin 
waves  are  present,  the  model  shows  the  approximate  additivity  of  spin-wave 
energies  discussed  in  Section  III. 

In  addition,  certain  other  energy  levels  in  our  model  can  be  coiqpared 

t 

with  exact  calculations  for  the  corresponding  states  in  class  C.  For  this 
purpose  it  is  convenient  to  introduce  the  abbreviation 

p,  =  s/n 


(26.) 
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We  shall  he  interested  in  states  near  the  antiferromagnetic  ground  state; 
i.e.,  those  for  which  p1  is  small. 

Consider  the  state  with  the  lowest  energy  for  a  given  p^:  the 
pockets  are  empty  for  0  s  x  s  np^  and  n(l-p^)  s  x  s  tt,  whereas  in  the 
interval  np.^  s  x  <:  tt(I-p^)  every  other  pocket  contains  a  ball.  (Compare 
this  with  CELI,  (l6)).  Our  model  predicts  an  energy: 

6  *  IN  7  [i-  cos  ('jr f>, ) ]  +  E0  (27) 

or,  for  small  values  of  p  , 

e-E0  ~  N  T  (28) 

The  exact  energy  for  the  corresponding  state  in  class  C  (see  CELI,  Pig.  l) 
lies  below  the  prediction  (27)  for  finite  p^;  the  asymptotic  formula  (28) 
is,  however,  correct,  at  least  to  within  .5  per  cent.  See  CELI,  Eq.  (Jl).1^ 
For  a  given  p^,  the  state  with  maximum  energy  in  the  model  corresponds 
to  enqpty  pockets  in  the  interval  n(^-p1)  s  x  <;  n^p^),  whereas  in  the  rest 
of  the  interval  every  other  pocket  contains  a  ball.  (Compare  this  with 
CELI  (18).)  The  model  predicts  an  energy: 

E  =■  IN  T  s  •'  n  (tT  ^, )  +■  Et  (29) 

and  for  small  p^: 

E-£,  ~  (30) 

Again,  the  exact  energy  for  the  corresponding  state  in  class  C  (see  CELI, 

Pig.  1)  lies  below  the  prediction  (29);  but  the  asymptotic  formula  (30) 
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is  correct.  See  CELI,  Eq.  (58). 

A  "uniform  density"  state  shall  he  one  in  which  for  every  interval 
b  s  x  s  d  such  that  d-b»l/N,  a  total  of  (d-b)pN/rr  pockets  contain 
balls.  The  "density"  p  is  equal  to  £  for  the  ground  state  (in  which 
every  other  pocket  is  filled)  and  the  quantity  is  equal  to  ^-p.  For 
such  states  our  model  predicts  an  energy: 

£  =  4-A/  T +  E0  (31) 

The  exact  energy  for  the  corresponding  state  in  class  C  is  easily  obtained 
by  a  straightforward  modification  of  the  Hulthdn  integral  equation,  (ill. 35) 
in  reference  2,  with  the  result: 

E-E.  =  toj  (l+lp,)  (32) 

Clearly  (31)  is  larger  than  (32)  for  finite  p^,  and  asymptotically  equal 
to  (32)  as  p^  goes  to  zero. 

On  the  basis  of  the  three  categories  of  states  discussed  above — 
minimum  energy,  maximum  energy,  and  "uniform  density" --it  is  tempting  to 
speculate  that  for  any  given  state  the  balls-in-pockets  model  always  yields 
an  energy  higher  than  the  exact  energy.  Of  this  we  have  no  proof.  In  any 
case,  the  fact  that  the  model  gives  the  correct  asymptotic  form  for  the 
energy  when  p^  is  small,  for  ail  three  categories,  as  well  as  for  the  des 
Cloizeaux  and  Pearson  spin  waves,  gives  one  confidence  that  it  provides  a 
good  description  of  the  energy  level  structure  of  class  C  near  the  anti- 
ferromagnetic  ground  state. 
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VI.  Partition  Function  for  Class  C:  Balls-in-Pockets  Model 

The  balls -in-pockets  model  discussed  in  Section  V  appears  to  provide 
u  good  approximation  to  the  low-lying  energy  levels  of  class  C.  In  this 
section  we  3hall  use  the  model  to  estimate  the  contribution  of  the  levels 
in  class  C  to  the  total  partition  function  of  the  antiferromagnetic  chain. 
The  problem  becomes  that  of  computing  the  partition  function  for  a  system 
of  balls  in  pockets,  as  shown  in  Fig.  4,  subject  to  rules  1  and  2  of  Section 
V  and  with  a  potential  energy  given  by  Eq.  (25). 

One  approach,  motivated  by  a  procedure  used  with  the  Ising  linear 
chain, l!+  is  as  follows.  Let  the  pockets  be  numbered  consecutively  from  1 
to  N-l  and  denote  the  state  of  the  i'th  pocket  by  an  Index: 


'  0  if  the  pocket  i a  empty. 

1  if  the  pocket  contains  a  ball. 


The  energy  of  a  state  shall  be  set  equal  to: 


E  =  Z  Pi*.)  +  Pv-,  V(>V,) 


where 

E  (  Pi,  Pic)  =  Pi  V(.xJ  -b  f(Plj  pltl) 

=  -f  (o,/)  =  -f  £i ,o)  =  0 

•f  0;  l)  =  -t  CO 

*£  =  7t  l  /N  (34) 


A  state  for  which  two  adjacent  pockets  both  contain  balls  will  have  an 
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Infinite  energy  due  to  the  function  f  and  therefore  will  not  contribute 
to  the  partition  function.  The  partition  function  for  the  model,  Zl.,  is 
given  by: 


=  Z  Z  *  Z  ^pf-pl  I.  £(TijPiH)  +  pu  V(*lw)]}(35) 

n«o  PwH  ‘s/  J 

where  g  =  l/kT. 

Let  be  a  2  x  2  matrix  with  rows  and  columns  labeled  by  the  indices 

and  P1+1>  respectively: 


Q(il  = 


Mi  0 


where 


Mi  =  tX/>[-(3V(x0]  . 


Let  T  be  the  column  vector 


r  » 


Zip  [-(3V(^.,)]  /  . 


Eq.  (35)  may  now  be  rewritten  as: 


-  Z  <?"'•  <? 


n*‘>  -  . 
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where  the  products  are  to  be  understood  as  matrix  products  and  E  denotes 
the  sum  of  the  two  elements  in  the  resulting  column  vector. 

If  the  matrices  were  all  equal  to  the  same  matrix  Q  with 

eigenvalues  of  unequal  magnitude,  the  evaluation  of  (39)  would  be  straight¬ 
forward.  One  could  diagonalize  Q  by  means  of  a  similarity  transforma¬ 
tion: 

q-  sv  S'1 


I a|  >  M 


The  product  of  the 


would  be  Q  raised  to  the  (N-2)nd  power: 


(40) 


1  c  tvV-2.  c  "I 

Q  -  SD  5  (41) 

But  for  large  N,  | \  lN-2»H  ,  and  the  partition  function  would  be 
equal  to  the  (N-2)nd  power  of  the  largest  eigenvalue  of  Q,  times  a  constant 
independent  of  N. 

In  our  case,  unfortunately,  the  matrices  depend  on  the  index  1. 

Nevertheless,  the  are  nearly  constant  for  i  restricted  to  a  range 

m  s  i  s  n  where  n-m«  N.  If  N  is  very  large  we  can  also  have  n-m»  1. 

The  considerations  of  the  preceding  paragraph  then  make  it  reasonable  to 
assume  that  in  the  limit  of  large  N  the  partition  function  (39)  is  equal  to: 


H-l 

B  TT 


(te) 
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where  B  is  a  constant  independent  of  N  and  A*  is  the  larger  (in  magnitude) 
of  the  two  eigenvalues  of  . 

We  have  given  an  intuitive  justification  for  (42);  in  fact,  it  is  not 
difficult  to  prove  the  result  rigorously  under  certain  restrictions  (which 
the  satisfy)  on  the  class  of  matrices  involved.  Since  the  proof  is  a 

bit  detailed  and  adds  no  further  insight  into  the  problem,  we  shall  not 
reproduce  it  here. 

The  largest  eigenvalue  of  (36)  is: 

“  i  +  if*  +  */*»)  *  (43) 


whence 

Z  iL,+ 

i-i 


At  this  point  it  is  convenient  to  introduce  the  renormalized  partition 
function  Zm  which  differs  from  in  that  all  energies  sure  measured  from 
the  ground  state  of  the  model.  Thus  Z^  should  provide,  at  low  temperatures, 
a  good  approximation  to  Z^,  the  contribution  of  states  in  class  C  to  the 
partition  function  of  the  linear  chain  when  energies  are  measured  from  the 
exact  ground  state.  Using  (26)  we  obtain  immediately: 


-  W  fc  =  h  z« 


(45) 
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f  * 


which  becomes  for  low  temperatures  (large  g);1^ 


The  normalized  temperature  'J'  is  equal,  to  (gj)-1.  For  small  the  free 
energy  (per  spin)  corresponding  to  (46)  is: 

rx/to  W 

The  coefficient  of  7*in  (47)  is  significantly  smaller  (in  magnitude) 
than  the  coefficient  in  (24).  In  our  opinion  the  disagreement  indicates 
that  the  overcounting  of  states  in  the  spin-wave  model  used  in  Section  IV 
leads  to  a  significant  error. 

The  quantity  Fg,  determined  by  numerical  integration  of  (44),  is 
plotted  in  Fig.  5  together  with  the  free  energies  for  short  chains  con¬ 
taining  10  and  11  spins.  This  figure  together  with  the  observation  that 

o 

the  coefficient  of  'f  in  (47)  is  only  half  the  coefficient  in  (7)  lends 
support  to  the  conclusion  that  even  at  low  temperatures  the  actual  free 
energy  of  the  antiferromagnetic  linear  chain,  F,  is  strictly  less  than 
(in  absolute  magnitude  strictly  greater  than)  F^.  If  this  last  inference 
is  correct,  the  discussion  at  the  beginning  of  Section  IV  allows  us  to 
conclude  that  makes  a  negligible  contribution  to  the  partition  function 
as  N  becomes  infinite.  Further  evidence  for  this  result  is  found  in 


N  _  f  *°  y  __  N 


Section  VII. 
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VII.  Class  C  and  the  Magnetic  Susceptibility 

In  CELI,  Section  IV,  ve  computed  the  magnetic  susceptibility  of  the 
antiferromagnetic  chain  at  T=0  utilizing  only  states  in  class  C,  and  only 
a  small  fraction  of  such  states.  The  method  there  employed  cannot  be 
extended  to  finite  temperatures.  In  fact,  by  considering  the  contribution 
of  the  states  in  class  C  to  the  susceptibility  at  finite  temperatures,  one 
obtains  rather  conclusive  evidence  that  Zc  is  a  negligible  fraction  of  the 
total  partition  function  for  long  chains. 

g 

The  number  of  states  in  class  C  with  total  spin  S  is : 


n  (S)  = 


(48) 


including  the  degeneracy  factor  2S+1.  Define: 

<s>c  =  (Zce-(!E"  0/(2* '■‘“"J  (.,) 

where  Sq  is  the  total  spin  and  Eq  the  energy  of  the  n'th  state,  and  the 
summation  goes  over  all  states  in  class  C.  The  fact  that  n(s)  is  very 
sharply  peaked  around  S=N/(2/5)  leads  one  to  suspect  that 

<S>c°c  N  (50) 

On  the  other  hand,  the  magnetic  susceptibility  (per  spin)  of  the 
chain  in  the  limit  of  zero  magnetic  field  is  given  by 

X  -  A  <  $X>/N 


(51) 
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where  A  is  a  constant  independent  of  N,  and  the  angular  brackets  denote 
an  average  as  in  (49),  but  with  summation  extending  over  all  energy  levels 
of  the  chain.  If  (S  )  in  (51)  were  replaced  by  (S  )c,  one  would  expect, 
on  the  basis  of  (50),  to  findX°CN;  i.e.,  an  infinite  susceptibility  for  a 
chain  of  infinite  length. 

The  foregoing  considerations  may  be  made  more  precise  as  follows. 

The  calculations  in  CELI,  Section  III  (compare  with  Section  V  above)  indi¬ 
cate  that  the  energies  of  all  states  in  class  C  with  total  spin  S  lie 
within  the  bounds 

0*  E  £  l  ~  TS  (52) 


For  the  considerations  which  follow,  the  weaker  condition 

0  <  E  i  M  f  ( S/N )  (53) 


is  sufficient,  where  f  denotes  any  positive,  continuous,  monotone  func¬ 
tion  of  x  such  that  f(0)=O. 

Let  us  assume  that  in  some  temperature  region  there  is  a  constant  p, 
(independent  of  N)  such  that 


Zc  /  H  i  X  >  0 


(5M 


where  Z ^  is  the  contribution  of  states  in  class  C  to  the  total  partition 
function.  Further  assume  that  the  susceptibility  is  finite  in  the  same 
temperature  region  where  (54)  holds: 

X  <  A, 


(55) 
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where  A^  is  a  constant  independent  of  N.  We  shall  now  show  that  the 
inequalities  (53)*  (5*0,  an4  (55)  lead  to  the  result  that  the  free  energy 
of  the  chain  is  constant,  and  therefore  the  entropy  is  equal  to  zero,  in 
the  temperature  region  for  which  (54)  and  (55)  hold.  Since  this  result 
is  highly  unlikely,  we  conclude  that  one  of  the  three  inequalities  from 
which  it  is  derived  must  be  invalid.  Of  these  three,  the  most  likely  to 
be  at  fault  is  (54).  That  is,  it  appears  highly  probable  that  the  contri¬ 
bution  of  states  in  class  C  to  the  total  partition  function  is  negligible 
as  N  becomes  infinite. 

The  proof  of  the  result  claimed  in  the  preceding  paragraph  is  obtained 
as  follows.  Since  both  the  quantities 

ZcS»e'<JE”  w  Z,S.e-fr* 


(where  the  second  sum  extends  over  all  states  in  class  B)  sure  positive,  the 
inequality  (54)  inqalies  that 

<S>  >  M  <S>C  /  C / +»)  (56) 

From  (51),  (55),  (56)  and  the  fact  that  (S2)  z  <S>2  it  follows  that 


<S>C  <  x  (57) 

where  is  independent  of  N.  Define 

T  -  XAxN^  (58) 

and  let  2L  and  7*  denote  sums  over  states  for  which  S  is  less  than  T  or 

C  C 


27 


greater  than  T,  respectively. 

A  consequence  of  these  definitions  and  Eq.  (49)  is  the  ineqjoality: 


(59) 


Then  (57),  (58),  and  (59)  imply  that 


Zc  e 


-|S£, 


"  >  Tc 


(60) 


or 


I',  c'**’  ^  lc  i  1  l'c 


(61) 


Define  n  by: 
o 


",  = 


I  H(S) 

5  =  / 


(62) 


Combining  (53),  (54),  and  (6l)  one  obtains  (note  that  has  been  replaced 
by  Z): 


-i*f  (*)  +  i  h "•  -  ii  z  i  £  h  £  +  v  h  "• 


(63) 


We  may  make  a  crude  estimate: 


I  £  ",  ^  N  nCr) 


(64) 


»  • 
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By  use  of  Stirling's  approximation  for  n(T)  and  the  ineqjoalities  (64), 
one  may  easily  verify  that  (l/N)  log  nQ  goes  to  zero  as  N  becomes  infinite. 
The  inequalities  (63)  then  imply  that  the  free  energy  is  identically  zero 
in  the  limit  of  infinite  N.  This  completes  the  proof. 
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VIII.  Conclusion 

The  balls -in-pockets  model  introduced  in  Section  IV  appears  to  provide 
a  very  good  description  of  the  energy  level  structure  in  class  C  near  the 
antiferromagnetic  ground  state.  Unfortunately,  this  does  not  enable  one 
to  conclude  very  much  about  the  low  temperature  thermal  properties  of  the 
chain,  because  the  partition  function  is  dominated  by  the  contributions 
from  levels  in  class  B.  The  upper  bound  on  the  free  energy,  (47),  may  not 
be  entirely  useless.  For  example,  it  rules  out  the  possibility  that 

-F  cc  TlA  (65) 

at  very  low  temperatures,  a  form  suggested  by  Eq.  (6). 

Equations  (47)  and  (7),  taken  together,  indicate  that 

lim  Ui  Zc  v  //* 

o  N  N ~¥  N  (66) 


for  T  >  0.  The  calculations  in  Section  VII,  on  the  other  hand,  show  that 


I  *•**  ?£  *  Q 

N-*  00  Z 

for  T  >  0,  on  the  basis  of  very  plausible  assumptions.  Strictly  speaking, 
the  result  (66)  is  stronger  than  (67).  But  the  arguments  leading  to  (66) 
are  much  more  tenuous;  in  particular  (7)  is  based  on  an  extrapolation  of 
results  for  short  chains  at  temperatures  on  the  order  of  |j/k  and  above. 

In  any  case,  the  conclusion  is  the  same:  namely,  an  adequate  discussion 
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of  the  thermal  properties  of  the  anti ferromagnetic  chain  at  low  tempera¬ 
tures,  if  based  upon  the  Be the  formalism,  must  take  into  account  the 
"bound"  states  in  class  B. 
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Appendix 

In  Sections  V  and  VI  we  disregarded  a  certain  difference  in  the 
counting  between  class  C  and  the  "balls -in-pockets"  model.  Associated 
with  a  single  state  in  the  model  is  a  set  of  2S+1  states  in  class  C 
forming  a  single  multiplet  with  total  spin  S.  It  is  easy  to  show  that 
this  "counting  error"  has  no  effect  upon  the  thermal  properties  in  zero 
magnetic  field. 

Let  Zg  be  the  contribution  of  states  with  a  given  total  spin  S  to 

the  partition  function  of  the  model.  Thus  the  model  partition  function 

Z  is  given  by: 
n 

W/a. 

Z-  -  Z  Z5  (• 

S*0 


If,  on  the  other  hand,  we  used  the  "correct"  counting  of  levels,  the 
partition  function  would  be 


z'  = 

»h- 

z  Us  -hi)  £s 

(A2) 

S-0 

Clearly, 

Zm  * 

Z'„  i  (.v+1)  Z „ 

(A3) 

and  hence 

h » 

■7  / 

_  Urn  ,09 

N-*o»  j\J 

(A4) 

That  is,  in  the  limit  of  large  N  the  free  energy  corresponding  to  Z^  is 
identical  with  that  corresponding  to  Zffl. 


Figure  Captions 


Fig.  1.  Energy  vs.  entropy  for  the  antiferromagnetic  chain.  The 
solid  curves  are  from  calculations  on  finite  chains;  the  dotted  curve 
represents  Eq.  (5).  The  dot-dash  curve  is  the  prediction  of  spin-wave 
theory  (reference  9)* 

Fig.  2.  Energy  vs.  temperature  for  the  antiferromagnetic  chain. 

The  solid  curves  are  from  calculations  on  finite  chains;  the  dotted  curve 
represents  Eq.  (6).  The  dot-dash  curve  is  the  spin-wave  prediction 
(reference  9). 

Fig.  3.  Some  sets  of  integers  for  the  case  N=l6  (see  text). 

Fig.  4.  Balls -in-pockets  model  for  N=l6  (see  text). 

Fig.  5-  Free  energy  vs.  temperature  for  the  antiferromagnetic  chain. 
The  dotted  curve  is  Fc,  determined  by  Eq.  (44). 
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